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Abstract 

We establish circumstances under which the dispersion of passive contami¬ 
nants in a forced, deterministic or random, flow can be consistently interpreted 
as a Markovian diffusion process. In case of conservative forcing the repul¬ 
sive case only, F = VV with V (x, t) bounded from below, is unquestionably 
admitted by the compatibility conditions. A class of diffusion processes is 
exemplified, such that the attractive forcing is allowed as well, due to an ap¬ 
propriate compensation coming from the ’’pressure” term. The compressible 
Euler flows form their subclass, when regarded as stochastic processes. 


Whenever one tries to analyze random perturbations that are either superim¬ 
posed upon or intrinsic to a driving deterministic motion, quite typically a con¬ 
figuration space equation x = v(x, t ) is invoked, which is next replaced by a for¬ 
mal infinitesimal representation of an Ito diffusion process dX(t) = b(X(t),t)dt + 
V2DdW(t). Here, W(t) stands for the normalised Wiener noise, and D for a diffu¬ 
sion constant. 

The dynamical meaning of b(x, t ) relies on a specific diffusion input and its pos¬ 
sible phase-space (e.g. Langevin) implementation, that entail a detailed functional 
relationship of v(x, t) and b(x. t), and justify such notions like: diffusion in an ex¬ 
ternal force held, diffusion along, against or across the deterministic how, [|l| , also 
with shear effects, ||. 
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The pertinent mathematical formalism corroborates both the Brownian motion 
of a single particle and the diffusive transport of neutrally buoyant components in 
flows of the hydrodynamic type. 

Clearly, in random media that are statistically at rest, diffusion of single tracers 
or dispersion of pollutants are well described by the Fickian outcome of the molecular 
agitation, also in the presence of external force fields (then in terms of Smoluchowski 
diffusions). On the other hand, it is of fundamental importance to understand how 
statistically relevant flows in a random medium (fluid, as example) affect dispersion. 
In the context of fluids, we might refer to diffusion enhancement due to turbulence, 
behaviour of Brownian particles in shear flows, but also to general effects of the 
external forcing (various forms of deterministic or random ’’stirring” of the random 
medium) exerted upon gradient or non-gradient, compressible and incompressible 
flows, and carried by them passive constituents, 

Except for suitable continuity and growth restrictions, necessary to guarrantee 
the existence of the process X(t) governed by the Ito stochastic differential equation, 
the choice of the driving velocity field v(x, t) and hence of the related drift b(x, t ) is 
normally regarded to be arbitrary. 

However, the situation looks otherwise, J2|. if we are interested in a diffusion 
of passive tracers in the a priori given flow whose velocity field is a solution of the 
nonlinear partial differential equation, be it Euler, Navier-Stokes, Burgers or the like. 
An implicit assumption, that passively buoyant in a fluid tracers have a negligible 
effect on the flow, looks acceptable (basically, in case when the concentration of a 
passive component in a flow is small). Then, one is tempted to view directly the fluid 
velocity field v(x, t ) as the forward drift b(x, t ) of the process, with the contaminant 
being diffusively dispersed along the streamlines. 

However, in general, the assumed nonlinear evolution rule for v(x, t ) must be 
checked against the dynamics that is allowed to govern the space-time dependence 
of the forward drift field b(x, t), ||, which is not at all arbitrary. The latter is ruled 
by standard consistency conditions that are respected by any Markovian diffusion 
process, and additionally by the rules of the forward and backward Ito calculus, 
0 , 1 - 

This particular issue we have analyzed before in the context of Burgers flows, 
||, where the Burgers velocity field was found to be inappropriate to stand for the 
forward drift of a Markovian diffusion process. Actually, the backward drift was a 
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correct identification. Then, the forced Burgers dynamics 

d t v B + (vb ■ V)Tb = DAv b + Vfi 


( 1 ) 


and the diffusion-convection equation 

d t c + [vb ■ V)c = DAc (2) 


for the concentration c(x, t ) of a passive component in a flow, in case of gradient 
velocity fields, were proved to be compatible with the Markovian diffusion process 
input. 

According to Ref. [J4|, in that case the dynamics of concentration results from 
the stochastic diffusion process whose density p(x, t ) evolves according to 


d t p = -DAp - V ■ (v B p) , 


(3) 


or equivalently: 

d t p=DAp-V -(bp) , (4) 

b = vb + 2 DVlnp . 

The previous reasoning can be easily exemplified by considering the standard un¬ 
forced Brownian motion with the initial (arbitrary, but sufficiently regular) density 
po(x). Its evolution po(x) —>■ p(x,t ) is implemented by the conventional heat kernel 
p(y,s,x,t ) = [4-7T D(t — g)]" 1//2 ea? 2 J>[— ]• The backward drift of the process (a 

solution of the unforced Burgers equation, originally denoted b*(x,t) in Ref. [||]) is 
defined as follows: Vb(x, t) = —2 DVlnp. The pertinent concentration dynamics is 
given by 

c(x,t) = f P*{y,0,x,t)c o (y)d 3 y (5) 

P*(y, 0, X, t) = p(y, 0, X, t ) . 

P{x,t) 

The remaining part is to determine the function Co(x,t ) i.e. the concentration 
of a tagged population in a Brownian ensemble. If we arbitrarily decompose the 
density of the process into p = pi + P 2 and regard pi(x, t ) as the density of a tagged 
population, then an appropriate definition of the concentration comes through: 


c(x, t ) 


Pi(x, t ) 

p(x, t) 


( 6 ) 
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By inspection one can check the validity of the diffusion-convection equation for 
c(x,t ) in a Brownian flow with the (backward drift) velocity v B (x, t) = —2 DWlnp. 

By combining intuitions which underly the self-diffusion description, with 
those appropriate for probabilistic solutions of the so-called Schrodinger boundary- 
data and next-interpolation problem, [£|. 0], the above argument can be generalized 

to conservatively forced diffusion processes. 

Namely, let us consider again the density p(x,t),t > 0 of a stochastic diffusion 
process, solving the Fokker-Planck equation d t p = DAp — V ■ (bp), where b(x,t ) 
stands for a forward drift. In case of conservative forcing, the drift solves an evolution 
equation: 

d t b + (b • V)6 = -DAb + VO . (7) 

For drifts that are gradient fields, the potential hi, whatever its functional form is, 
must allow for a representation formula, reminiscent of the probabilistic Cameron- 
Martin-Girsanov transformation: 

(l(x,t) = mdfi + l(A + V-b)} , (8) 

where b(x,t ) = V3>(x, t). 

For the existence of the Markovian diffusion process with the forward drift b(x, t ), 
we must resort to potentials Q(x,t) that are not completely arbitrary functions. 
Technically, jSJ, the minimal requirement is that the potential is bounded from 
below. This restriction will have profound consequences for our further discussion 
of diffusion in a flow. 

If we set p = pi + p 2 again, and demand that pi ^ p solves the Fokker-Planck 
equation with the very same drift b(x, t ) as p does, then as a necessary consequence 
of the general formalism, [§§• the concentration c(x, t) = solves an associated 

diffusion-convection equation d t c+ (v B ■ V)c = DAc. Here, the flow velocity v B (x, t ) 
coincides with the backward drift K=v B of the generic diffusion process with the 
density p(x,t ) and reads: v B = b — 2D'Vlnp. Obviously, the forced Burgers equation 
(1) is identically satisfied. 

We should clearly discriminate between forces whose effect is a ’’stirring” of the 
random medium and those acting selectively on diffusing particles, with a negligible 
effect on the medium itself. For example, the traditional Smoluchowski diffusion 
processes in conservative force holds are considered in random media that are sta¬ 
tistically at rest. Following the standard (phase-space, Langevin) methodology, let 
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us set b(x) = jjK(x), where j3 is a (large) friction coefficient and K represents an 
external Newtonian force per unit of mass ( e.g. an acceleration) that is of gradient 
from, K = —VC/. Then, the effective potential V reads: 


Q 


K 2 
2/P 




(9) 


and the only distinction between the attractive or repulsive cases can be read out 
from the term V-/?. For example, the harmonic attraction/repulsion K = =|=q:x, a > 
0 would give rise to a harmonic repulsion, if interpreted in terms of Vfi, in view of 
Q = =F 3D^. The innocent looking =)=3 D^j renormalisation of the quadratic 

function gives rise to entirely different diffusion processes, with an equilibrium mea¬ 
sure arising in case of U(x) = +^x 2 only. 

The situation would not change under the incompressibility condition (cf. also 
the probabilistic approaches to the Euler, Navier-Stokes and Boltzmann equations, 
0). Following Townsend’s, [|2|, early investigation of the diffusion of heat spots in 
isotropic turbulence we may choose U(x) = ^x 2 — j(y 2 + z 2 ) which implies V- K = 0. 
Then, V(x) = f^[x 2 + \(y 2 + £ 2 )], hence the repulsive is produced again in the 
equation of motion characterising a stationary diffusion in an incompressible fluid: 
div v = 0, b = b* = v —> (v ■ V)F = Vfh 

By formally changing a sign of we would arrive at the attractive variant of the 
problem, that is however incompatible with the diffusion process scenario in view of 
the unboundedness of —V from below. 

We have thus arrived at the major point of our discussion: we may get in trou¬ 
ble with the Markovian diffusion input in case of general external ’’stirring” forces. 
Hence, we must specify an admissible class of perturbations which, while modify¬ 
ing the flow dynamics, would nonetheless generate a consistent diffusion-in-a-flow 
transport of passive tracers. 

Should we a priori exclude the attractive variants of the potential ? Can 
we save the situation by incorporating, hitherto not considered, ’’pressure” term 
effects as suggested by the general form of the compressible Euler (here F = — W 
stands for external volume forces and p for the fluid density that itself undergoes a 
stochastic diffusion process): 


d t v E + (ve ■ V)Fb = F - -VP 

P 


( 10 ) 
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( 11 ) 


or the incompressible, 0, Navier-Stokes equation: 

dtVNs + (vns ■ V)vns = —Avns + F -VP , 

P P 

both to be compared with the equations (1) and (7), that set dynamical constraints 
for respectively backward and forward drifts of a Markovian diffusion process ? 

Notice that the acceleration term F in equations (10) and (11) normally is re¬ 
garded as arbitrary, while the corresponding term Vf2 in (1) and (7) involves a 
bounded from below function f2(x, t). 

Since, in case of gradient velocity holds, the dissipation term in the incompress¬ 
ible Navier-Stokes equation (11) identically vanishes, we should concentrate on an¬ 
alyzing the possible ” forward drift of the Markovian process” meaning of the Euler 
how with the velocity held ve, (10). 

At this point it is useful, at least on the formal grouds, to invoke the standard 
phase-space argument that is valid for a Markovian diffusion process taking place 
in a given how v(x, t ) with as yet unspecihed dynamics. We account for an explicit 
force exerted upon diffusing particles, while not necessarily directly affecting the 
driving how itself. Namely, [j|j. S], let us set for inhnitesimal increments of phase 
space random variables: 

dX(t) = V(t)dt 

dV(t) = (3[v{x, t) - V(t)]dt + K{x)dt + (dV2DdW(t) . (12) 

Following the leading idea of the Smoluchowski approximation, we assume that 
(3 is large, and consider the process for times significantly exceeding /3 _1 . Then, an 
appropriate choice of the velocity held v(x, t ) (boundedness and growth restrictions 
are involved) may in principle guarrantee, [;3l], the convergence of the spatial part 
X(t) of the process (12) to the Ito diffusion process with inhnitesimal increments: 

dX(t) = v(x, t)dt + V2 DdW(t) . (13) 

However, one cannot blindly insert in the place of the forward drift v(x, t) any 
of the previously considered bulk velocity helds, without going into apparent con¬ 
tradictions. Specifically, the equation (7) with v(x,t) b(x,t ) must be valid. 

By resorting to velocity helds v(x, t) which obey Av(x, t) = 0, we may pass from 
(7) to an equation of the Euler form, (10), provided (8) holds true and then the 
right-hand-side of (7) involves a bounded from below effective potential fh 
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An additional requirement is that 


F — -VP = vn . 
P 


( 14 ) 


(b 


Clearly, in case of a constant pressure we are left with the dynamical constraint 
^ ve)' 

d t b + (b- V)b = F = Vft (15) 


combining simultaneously the Eulerian fluid and the Markov diffusion process in¬ 
puts, if and only if F is repulsive, e.g. — V (x, t ) is bounded from below. Quite 
analogously, by setting F = 0, we would get a constraint on the admissible pressure 
term, in view of: 

d t b+ (b- V)6 = —-VP = VO . (16) 

P 

Both, in cases (15), (16) the effective potential O must respect the functional 
dependence (on a forward drift and its potential) prescription (8). In addition, the 
Fokker-Planck equation (4) with the forward drift Ve{x. t) = b(x,t ) must be valid 
for the density p(x, t ). 

To our knowledge, in the literature there is known only one specific class of 
Markovian diffusion processes that would render the right-hand-side of Eq. (10) 
repulsive but nevertheless account for the troublesome Newtonian accelerations, e.g. 
those of the from — W, with +V bounded from below. Such processes have forward 
drifts that for each suitable, bounded from below function V(x) solve the nonlinear 
partial differential equation: 


d t b+(b-V)b=-DAb + V(2Q-V) (17) 

with the compensating pressure term: 

Q = 2-D 2 —T 7 y- = \u 2 + DV • u (18) 

p 2 

u(x, t ) = DWln p(x, t ) 

Their exhaustive discussion can be found in Refs. BIBB, together with indi¬ 
cations for their possible relevance as a stochastic counterpart of the Schrodinger 
picture quantum dynamics. Clearly, we have: 

F = -VE , V2Q = --VP (19) 

P 
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where: 


P(x, t ) = —2 D 2 p(x, t ) A In p(x, t ) 


( 20 ) 


Effectively, P is here dehned up to a time-dependent constant. Another admissible 
form of the pressure term reads (summation convention is implicit): 


A;. [p ( 2DV,.t }h: jt\ = V,(2Q) 


( 21 ) 


If we consider a subclass of processes for which the dissipation term identically 
vanishes ( a number of examples can be found in Refs. [S|): 


A b(x, t) = 0 


( 22 ) 


the equation (17) takes a conspicuous Euler form (10), ve b. 

Let us notice that (20), (21) provide for a generalisation of the more familiar, 
thermodynamically motivated and suited for ideal gases and fluids, equation of state 
P ~ p. In case of density fields for which — Ain p ~ const , the standard relationship 
between the pressure and the density is reproduced. In case of density fields obeying 
—A In p = 0, we are left with at most purely time dependent or a constant pressure. 
Pressure profiles may be highly complex for arbitrarily chosen initial density and/or 
the flow velocity fields. 

To conclude the present discussion let us invoke Refs. SSI- The problem of 
a diffusion process interpretation of various partial differential equations has been 
extended beyond the original parabolic equations setting, to nonlinear velocity held 
equations like the Burgers one, see e.g. [3|]. On the other hand, the nonlinear 
Markov processes associated with the Boltzmann equation, in the hydrodynamic 
limit, are known to imply either an ordinary differential equation with the velocity 
held solving the Euler equation, or a diffusion process whose drift is a solution of 
the incompressible Navier-Stokes equation (without the curlv = 0 restriction), [^]. [8|. 
The case of external forcing has never been satisfactorily solved. 

Our reasoning went otherwise. We asked for the admissible space-time depen¬ 
dence of general velocity helds that are to play the role of forward drifts of Markovian 
diffusion processes. Our finding is that solutions of the compressible Euler equation 
are appropriate for the description of a non-deterministic (e.g. random and Marko¬ 
vian) evolution and belong to a class of Markovian diffusion processes orginally 
introduced by E. Nelson in his quest for a probabilistic counterpart of the quantum 


8 


dynamics, [p]. [6J. Our solution of the problem involves only the gradient velocity 
fields. However, a couple of issues concerning the curl b ^ 0 velocity fields and their 
nonconservative forcing have been raised in Refs. [[J. 
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